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We discuss the spatial structure of the Cooper pair in dilute neutron mat- 
ter and neutron-rich nuclei by means of the BCS theory and the Skyrme- 
Hartree-Fock-Bogioliubov model, respectively. The neutron pairing in 
dilute neutron matter is close to the region of the BCS-BEC crossover in 
a wide density range, giving rise to spatially compact Cooper pair whose 
size is smaller than the average interaparticle distance. This behavior ex- 
tends to moderate low density (~ 10" 1 of the saturation density) where 
the Cooper pair size becomes smallerst (~ 5 fm). The Cooper pair in 
finite nuclei also exhibits the spatial correlation favoring the coupling 
of neutrons at small relative distances r^3 fm with large probability. 
Neutron-rich nuclei having small neutron separation energy may provide 
us opportunity to probe the spatial correlation since the neutron pairing 
and the spatial correlation persists also in an area of low-density neutron 
distribution extending from the surface to far outside the nucleus. 



1. Introduction 



The formation and the condensation of the Cooper pairs are the essence of 
superconductivity and superfluidity in many-Fermion systemsPJ The bind- 
ing energy of the Cooper pair is closely related to the pairing gap A. The 
spatial size of the Cooper pair is identified to the coherence length £ of 
the superconductors, which plays important roles in many aspects, for in- 
stance, in distinguishing the type I and type II superconductors. What is 
the size of the Cooper pair in the superfluidity of nuclear systems? A simple 
estimate of the coherence length £, based on the uncertainty principle in 
uniform matter, leads to £ ~ with Dp being the Fermi velocity^ If one 
considers saturated nuclear matter as a simplification of finite nuclei, and 
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adopts the typical value of the pairing gap A ps 12 /y/A ~ 1 MeV appropri- 
ate for heavy nuclei, the estimate gives £ ~ 20 fm which is much larger^! 
than the radius of nuclei R « 1.2 A 1 / 3 ~ 3 — 7 fm or interparticle distance 
~ 2.5 fm in saturated matter. However, if one considers extreme situations, 
such as dilute neutron matter and exotic nuclei with large neutron excess, 
there appear new features of the nuclear pairing that can be related to the 
spatial structure of the Cooper pair. It is the aim of this article to illustrate 
it using a few examples. 

2. Dilute neutron matter 

The superfluidity in neutron matter is density dependent.^^ The pairing 
gap can be obtained by solving the BCS equations for the bare nuclear 
force in the X S channel at each neutron density p = kp/3w 2 or the Fermi 
momentum kp. The gap is small A <C 1 MeV at hp = 1.36 fm -1 (p/po = 
1, the neutron density at saturation po = 0.08 fm -3 ). With decreasing 
the density it first increases, reaching the maximum A s» 3 MeV around 
kp pa 0.8 fm (p/po ~ 0.2), then decreases and approaches to zero at the 
low-density limit. Other many-body medium effects which are beyond the 
BCS approximation reduce the gap, but the predictions vary depending on 
the theoretical methods.®^ Recent ab initio Monte Carlo calculations E2HH 
on the other hand, predict rather modest reduction by less than 50%, and 
the qualitative features of the density dependence is kept. Having these 
reservations in mind, let us consider the structure of the neutron Cooper 
pair in the BCS approximation.^ 

The Cooper pair wave function can be defined, apart from the normal- 
ization, as an expectation value of the pair operator with respect to the 
BCS state: 

% air (ri, r 2 ) = (ip(r! tW(r 2 I)) = £ u k v k e ik - r . (1) 

k 

It is a function of the relative coordinate r = r 2 — i*i of the two neu- 
trons, and in the momentum space it is a product of the u and v factors. 
Examples of the Cooper pair wave functions are shown in Fig. 1 for two 
different densities.^ The wave function exhibits an oscillatory behavior 
characterized by the Fermi wave length 2n/kp and an overall decay profile 
whose asymptotic form is exponential ~ exp(— rA/hvp) ( for large relative 
distance r = |r 2 — r\\) whose length scale is nothing but the coherence 
length, or the size of the Cooper paper.E More precisely, the coherence 
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length can be calculated as the rms radius of the Cooper pair £ = 




with (r 2 ) = J drr 2 \V pair (r)\ 2 / J dr\V pair (r)\ 2 . 

An interesting feature of the neutron Cooper pair in superfluid neutron 
matter is that its size also varies significantly with changing the neutron 
density (See Fig. 1(a)). From a very large value £ = 46 fm at p/po = 1, 
the coherence length £ decreases sharply with decreasing the density. The 
coherence length takes the smallest values £ = 5 — 8 fm for a rather wide 
range of the density p/po = 0.2 — 10~ 2 , and it increases gradually with 
decreasing the density. 




P'Po r (fm) 

Fig. 1. (a) The coherence length £ and the average interparticle distance d in superfluid 
neutron matter, plotted as a function of the neutron density p/po (with po = 0.08 fm~ 3 ). 

(b) (c) The Cooper pair wave function r 2 \ty pa i r (r)\ 2 at densities p/po = 1 and 1/8. The 
bare force is used in the present BCS calculation!^ 

The Cooper pair wave function at densities where the coherence length 
is the smallest is very different from that of of the electron Cooper pair 
in the traditional metal superconductors. An example is shown in Fig.l 

(c) , which is for p/po = 1/8 where the coherence length £ = 4.9 fm is 
close to the minimum value. It is seen that the oscillatory behavior is 
strongly suppressed. The probability distribution is concentrated (~ 80%) 
at small relative distances within the first node r < ir/kp w 4.5 fm, and the 
probability at the second and third bumps is very small. This is because 
the size of the Cooper pair (£ = 4.9 fm) is almost equal to the position of 
the first node n/kp which is nothing but the average interparticle distance 
d = p- 1/3 n/kp. The size of the Cooper pair is "small" in this sense. 
This is quite contrasting to the metal superconducters where the Cooper 
pair size £ is thousands times larger than the average interparticle distance 
d. The situation of the "small" Cooper pair £ ^ d is seen in a wide interval of 
densities p/po = I0~ 4 — 10 _1 (Fig. 1(a)). The shape of the Cooper pair wave 
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function at these densities is similar to that of Fig. 1(c), and the probability 
is even more concentrated in the first bump although the absolute size 
is larger at very low densities p/po ~ 10~ 4 — 10~ 2 . It is noted that the 
Cooper pair at moderate low densities p/po m 10 _1 — 0.5 exhibits also the 
strong spatial correlation at small relative distances. The wave function at 
p/po = 0.5 is shown in Fig. 1(b). In this case the calculated coherence length 
£ = 11 fm is a few times larger than the average interparticle distance d = 
2.8 fm. Nevertheless the concentration of the probability within the relative 
distance r^3 fm (in the first bump) is signfinicant, and the probability in 
r < 3 fm reaches as large as ~ 50%. 

The situation of the small Cooper pair £/d^l is related to the so- 
called BCS-BEC crossover phenomenonp*"^ which has been discussed in- 
tensively in ultra-cold Fermi atom gas in a trapPSEH It is a phenomenon 
which can occur generally in any kind of many-Fermion superfuluid sys- 
tems by changing the strength of the interparticle attractive force or the 
density. In a situation of the weak interaction, which the original BCS 
theory has dealt with, the bound pair (the Cooper pair) can be formed 
only in the medium. However, if the interaction is as strong as to form 
a bound pair (a composite boson) even in the free space, the condensed 
phase is more close to a condensate of the composite bosons, i.e. the Bose- 
Einstein condensate (BEC). The BCS-BEC crossover is characterized by 
the ratio £/d of the coherence length and the average interparticle distance 
and the ratio A/ep of the pairing gap and the Fermi energy. The weak- 
coupling BCS and the BEC limits correspond to £/d ~> 1, A/ep <C 1 
and £/<i <§; 1, A/e F ^ 1, respectively while the region of the crossover 
may be related to 0.2££/d£l.2 and 0.2 <> A/e F & 1.3PE 3 At the mid- 
way of the crossover, called the unitarity limit, the interaction strength is 
on the threshold to form the isolated two-particle bound state, and the 
values are £/d = 0.36, A/ep — 0.69. In the BCS calculation discussed 
above small £/d ratio 0.7 — 1.2 and large A/ep ratio 0.2 — 0.4 is realized 
at p/po ~ 10~ 4 — 10 _1 . (Note that also in an ab initio calculation,^ the 
large gap ratio A/ep ~ 0.2 — 0.3 is obtained in approximately the same 
but slightly small density region.) We can regard dilute neutron matter in 
the wide low-density interval p/po = 10~ 4 — 10 _1 (or in slightly narrower 
interval) as being in the crossover region. We note here that the nuclear 
force in the X S channel has a large scattering length a = — 18 fm, indicat- 
ing that the interaction strength is very close to the threshold to form a 
two- neutron bound state. The small Cooper pair £/d ^ 1 at low densities 
originates from the nature of the nuclear force. 
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3. Cooper pair in neutron-rich nuclei 

Let us consider the spatial structure of the Cooper pair in finite nuclei. 

The spatial structure of the correlated two neutrons has been discussed 
intensively for two neutrons in the light two-neutron halo nuclei 11 Li and 
6 He in (inert or active) core plus two neutron models!^!®] A common 
prediction is that the valence halo neutrons exhibit a spatial correlation 
favoring the 'di-neutron' configuration with two neutrons coupled at small 
relative distances. The spatial correlation is also discussed in stable heavy 
nuclei with closed-shell core plus two neutrons, e.g. 206 > 210 Pb, by means of 
shell model approaches pSHHH One can generalize these findings by using the 
Hartree-Fock-Bogoliubov (HFB) method, which can be applied to a wide 
class of open shell nuclei including isotopes very close to the drip-line and 
also to non-uniform matter. 

Let us start defining the wave function of the Cooper pair in finite nuclei. 
It may be given by 

*p«rir(ri,r2) = (*A-2W(ri t)V(r 2 1)\*a) (2) 

using the pair correlated ground states &a and 3>a-2- This represents the 
probability amplitude of removing two neutrons (positioned at r± and r%) 
from the ground state $,4, and leaving the remaining system in the ground 
state &a-2- Provided that the ground state is described within the HFB 
framework, where the ground states with different nucleon numbers are rep- 
resented by a single HFB state $hfb, the definition Eq.Q can be replaced 
with the expectation value as in Eq.([T]). Then, since the HFB state is a 
generalized Slater determinant consisting of the Bogolviubov quasiparticlc 
states, this quantity is evaluated^ as a sum over all quasiparticle states i 

%air(ri,r 2 ) = ($HFB|^(ri t)V>(r 2 |)|$hfb) = £ ^Vi tVf } >2 I) 

i 

(3) 

using the first and the second components of the quasiparticle wave func- 
tion (f>i(ra) = ((p^\ra), ip^ 2 \rcr)). In the following we show the results of 
our HFB calculation, which adopts the Skryme functional and the density- 
dependent contact interaction as a phenomenological pairing force 121123 The 
parameter set of the pairing interaction is such that it reproduces the scat- 
tering length a = — 18 fm in the low-density limit, and reproduces the 
average pairing gap in known nuclei 021221 

An example calculated for 142 Sn is shown in Fig. 2. Here one neutron is 
fixed at the position slightly outside the nucleus r% — 7 fm and the prob- 
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Fig. 2. (Left) Cooper pair wave function |$po*r(*"l) fi)\ 2 / Pni^l) i n neutron-rich nu- 
cleus 142 Sn, plotted as a function of V2 on the xz plane while the coordinate T\ is fixed 
to (0, 0, 7) fm located slightly outside the surface. (Right) The same but plotted along 
the z-axis. Different curves are results obtained by putting cut-off's with respect to the 
orbital angular momentum I of the single-particle orbits. 

ability distribution | ^ pa ir (fi, r 2)\ 2 is plotted as a function of r^. It shows 
that the second neutron has a large probability (~ 50%) to be correlated at 
small relative distances |r*x — »"2 1 3 fm to the partner neutron. The spatial 
correlation seen here is generic in a sense that it is seen systematically in 
Ca, Ni, and Sn isotopes including both stable and neutron-rich nuclei.^ 
The strong spatial correlation is also seen in other HFB calculations which 
adopt the finite-range Gogny force as the effective pairing forceP^ffl 

I emphasize here that a large single-particle space is necessary in de- 
scribing the spatial correlation.^ In order to describe the correlation 
with the length scale D ~ 3 fm, the single-particle basis needs to cover 
a momentum range up to p m ax ~ h/D, which corresponds to a maximal 
energy 6 max 

~ Pmax/^ m ~ 80 MeV, or a maximal angular momentum 
Imax ~ Rpmax ~ 10ft (for the nuclear radius R ~ 5 fm). This is demon- 
strated in Fig.2(right), where the summation over the quasiparticle states 
i in Eq. ((3]) is truncated by introducing a cut-off with respect to the orbital 
angular momentum I. Single-particle orbits with large angular momentum 
up to Imax ~ 10 have sizable contributions. Note that in 142 Sn with N = 92 
the Fermi energy is around the orbit, and the maximal orbital angular 
momentum of the orbits occupied in the independent particle limit is Z = 5. 
The single-particle states with I = 5 — 10 lie high above the Fermi energy. 
If one uses the harmonic oscillator basis, it should include ~ 10 oscillator 
quanta. In fact, all the HFB calculation ;] 36 ! 39 ^ where the strong spatial 
correlation in the Cooper pair wave functions is demonstrated adopt such 
a large single-particle space. 

Equivalently, a small single-particle space is insufficient. If we restrict 
ourselves to a single-j shell (nlj), i.e., the sum in Eq.([3]) is restricted to the 
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magnetic substates of the orbit (nlj), we obtain the angular correlatio n 1 41 ! 42 ! 
Pi (Q 12) for small relative angles #12 Ss l/l, but the correlation with respect 
to the radial direction is not produced. Inclusion of all the orbits in one 
oscillator shell still has deficiency! 39 ! 40 ! The Cooper pair wave function in 
this case exhibits an artificial symmetry ^ pa ir{i' 1^2) — ±.^ pa ir{— r 2 ) 
because of the common single-particle parity, and the probability appears 
not only around r 2 ~ r±, but also around the mirror reflected position 
r 2 ~ -rx. 



2.0 1 1 0.20 2.0e-6n-r 




Fig. 3. (a) Pair potential A(r) in 142 Sn obtained in the Skyrme-HFB model^llMl (solid 
curve), and the neutron density pn(r) (dotted curve), (b) The Cooper pair wave function 
\y p air(ri,r 2 )\ 2 in 142 Sn and 120 Sn, with the first coordinate fixed at ri = (0,0,9) fm 
and plotted as a function of 7*2 along the z-axis. 



The spatial correlation of neutron Cooper pairs plays an important role 
if we consider neutron-rich nuclei with small neutron separation energy. 
Nuclei of this kind often accompany low-density distribution of neutrons, 
called skin or halo, extending from the nuclear surface toward the outside. 
Figure 3(a) is an example of the pair potential A(r) for the very neutron rich 
nucleus 142 Sn obtained in the same Skyrme-HFB calculation as in Fig. 2. 
The pair potential A(r) exhibits significant enhancement around r ~ 5 — 8 
fm, which is slightly outside the nuclear surface (the corresponding neutron 
density there is about 1/2-1/10 of the central density). The pair potential 
decreases rather slowly with moving outside the surface region, and it is 
about to diminish only at very large distances r <; 12 fm. It is much more 
extended than the neutron density. Furthermore the spatial correlation 
persists in this low density region as shown in Fig. 3(b). We note that the 
spatial correlation is present also in stable open shell nuclei,^ 39 ! 40 ! and it is 
enhanced around the nuclear surface. However the nucleons (and hence the 
Cooper pairs) do not penetrate far outside the surface in stable isotopes (cf 
bottom panel of Fig. 3(b)). The pair correlations in the dilute surrounding 
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is a unique feature of weakly bound nuclei. 

4. Probing the spatially correlated Cooper pair 
4.1. Soft modes 

If spatially correlated di-neutrons exist in nuclei, especially in the low- 
density skin/halo region, there may emerge new modes of excitation re- 
flecting the motion of di-neutron(s). This simple ide a! 22 * 24 * has been a focus 
of theoretical and experimental studies of the soft dipole excitation in two- 
neutron halo nuclei. Although the reality is not that simple, the core+n+n 
models ^ * 27 * 30 ! of n Li explain the observed large El strength of the soft 
dipole excitation^ in terms of the pairing and the spatial correlation of 
the valence halo neutrons. It is interesting to explore possibility of similar 
excitation modes in heavier mass neutron-rich nuclei, where more than two 
weakly bound neutrons contribute to the pair correlation. 

A useful scheme to describe excitation modes built on the pair cor- 
related ground state is the quasiparticle random phase approximation 
(QRPA). Let us take the formulation based on the same Skyrme-HFB 
model that is used for the description of the ground state l221iiliSI Having a 
QRPA excited state \n,LM), one can calculate the two-particle amplitude 
(n 1 LM\il^{r\ t)V^( r 2 4)|0 9S ) which tells us how two particles move in the 
excited state \n, LM) in reference to the ground state (of the N — 2 sys- 
tem). For simplicity let us look at the zero-range part at ri = T2 of the 
amplitude: 

Pr r (r) = (n, LM\^(r t)^(r I)|0 9S ) , (4) 

which is called the pair transition density. 

Figure 4 is an example of soft dipole excitation which suggests motion of 
the spatially correlated di- neutrons.^ The soft dipole excitation is seen here 
as a bump of the El strength which lies just above the neutron separation 
energies (Si n , S^n = 1-9, 2.4 MeV). In neutron-rich Ni isotopes beyond the 
N = 50 shell closure both of the one- and the two-neutron separation 
energies are calculated to be very low S n , S^n ~ 1 — 3 MeV. In such weakly 
bound nuclei, the low-lying dipole modes appear just above the separation 
energy since it is possible to excite a bound neutron to unbound orbits in 
the continuum, letting the neutron escape from the nucleus. If the pair 
correlation is taken into account, however, the mode is dominated by the 
pair motion rather than by a simple particle-hole (or independent two- 
quasiparticle excitation). Consequently the pair transition density P^ air {r) 
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Fig. 4. (a) The B(E1) strength function in neutron-rich nucleus 84 Ni, obtained with the 
Skyrme-HFB + continuum QRPA methodfS] The large strengths around E = 10 — 20 
MeV are the giant dipole resonance (GDR) while some amount of strength is distributed 
just above the one- and two-neutron separation threshold energies (the small arrows), 
(b) The neutron pair transition density pP a,r (r) of the soft dipole mode evaluated at 
E = 3 MeV (marked with the big arrow in (a)). Figures taken from RefPH 



has larger amplitude, especially for r > R SU rf, as seen in Fig. 4(b). It is 
not explicit in this figure whether the neutron pair in the excited state 
is spatially correlated, but we can infer it from the observation that a 
large number of orbital angular momenta I reaching more than 10h have 
significant and coherent contributions to the pair transition density. As we 
discussed above (cf. Fig. 2), large I implies a spatial correlation at small 
distances between the two neutrons. 

A similar mode of excitation having the character of di-neutron motion 
is predicted also in the octupole response in the same istopes >80 Ni beyond 
N = 50P It is a smooth distribution of neutron strength lying just above 
the threshold energy (like the soft dipole mode), and it coexists from the 
octupole surface vibrational mode of the isoscalar character seen in many 
of stable nuclei. 

In contrast to the light two-neutron halo nuclei, the presence of the 
spatial correlation does not influence strongly the El strength of soft 
dipole excitation in heavy neutron-rich nuclei such as 84 Ni. We need 
other probes which are directly connected to the pair transition density. 
Since the soft dipole excitation in 11 Li and in >80 Ni is located above the 
two-neutron separation energy, one can expect that momentum distribu- 
tion/correlation of two neutrons emitted from the soft mode may carry 
information on the spatial correlation of the neutron pair. Quantitative 
theoretical description of the two-neutron correlation is achieved only for 
the corc+n+n models^^ for n Li and 6 He, and experimental information 
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is very scarce so far! 48 * 49 ! It is possible to describe the two-neutron cor- 
relation in the continuum also in the framework of the QRPA since the 
information on the directions of two neutrons are contained in the pair 
transition density (n, LM\ifft(ri t)V' t ( T *2 4-)|0 gs ) especially in its asymp- 
totic form at \r\\, | t*2 | —> oo. It is an interesting future subject to study in 
heavier neutron-rich nuclei such as >80 Ni using the HFB+QRPA formalism. 

4.2. Two-neutron transfer 

The two-neutron transfer reactions such as (p,t) and (t,p) are known as a 
good probe to the pair correlation in the ground stateP^* 50 * 51 ! More pre- 
cisely it can be regarded as a probe of the Cooper pair wave function, espe- 
cially its behavior at small relative distances between the paired neutrons. 
Consider the (p,t) reaction populating the ground state of the neighboring 
N — 2 nucleus in the single-step DWBA and the zero-range approximation. 
Then the transition matrix elements involves the form factoi^^l 

F(R)= J dr (0 gSjN _ 2 \ij(R + r/2t)ip(R- r/2 l)\0 gs , N ) <f>{r) 

which is the convolution of the Cooper pair wave function ty pa i r (R + 
r/2, R — r/2) with the two-particle wave function <f>{r) in the triton. Not- 
ing the small radius of the triton ~ 2 fm, we see immediately the form 
factor picks up the correlation at small relative distances in the Cooper 
pair wave functions. It is then not a very bad approximation to uti- 
lize the Cooper pair wave function at zero relative distance r = 0, i.e. 
1 i pa ir(R, R) = (ip(R I)) = Ppair(R) as a substitute of the form 

factor assuming F(R) cx P pa i r {R). P pa i r (R) is nothing but the pair den- 
sity p(R) implemented automatically in the Skyrme-HFB model using the 
pairing force of the contact typeP^ 

An example^ of the calculated pair transition density P pa i r (R) is shown 
in Fig. 5(a) for Sn isotopes covering from stable isotopes to very neutron-rich 
150 Sn. It is seen that the radial dependence of P pa i r {R) suddenly changes 
at the AT = 82 shell closure (at 132 Sn). In neutron- rich isotopes beyond 
N = 82, the amplitude extends far outside the nuclear surface r > R sur f+3 
fm (> 9 fm). This happens because neutron single-particle orbits above the 
N = 82 shell gap are bound only weakly, and the weakly bound neutrons 
have density distributions extended far outside the nuclear surface. (The 
one-neutron separation energy is the order of ~ 2 — 3 MeV for A > 132, 
while it is more than 8 MeV in isotopes with A < 132.) Consequently 
both the pair potential A(r) and the Cooper pair wave function keep non- 
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Fig. 5. (a) Neutron pair transition density P pa i r (R) for the ground state transition, 
evaluated for even-even Sn isotopes for A = 120 — 130(thin solid curve), A = 134 — 140 
(dotted), and A = 142 - 150 (thick solid), (b) The two-neutron transfer strength B(P0) 
for the ground state transition (filled diamond), the strength for two- neutron addition 
transfer for the excited + state (small triangle), and the squared pairing gap A 2 (open 
circle) for the even-even Sn isotopes. The horizontal axis is the mass number A. 



negligible magnitude even far outside (See also Fig. 3). As seen in the figure 
the amplitude P pair (R) extends up to r ~ 12 fm for the isotopes A > 140. 

The above observation leads to an expectation that the (p,t) and (t,p) 
cross sections may be enhanced considerably as the neutron separation en- 
ergy becomes smallP^ An estimate of the isotopic trend, much simpler than 
the DWBA calculation, is shown in Fig. 5(b). Here is plotted the 'strength' 
which is defined by B(P0) = J dRP pa i r (R)\ . It is illuminating to com- 
pare it with the isotopic trends of the pairing gap A squared (A being an 
average value of the pair potential A(r)). If the pair potential and the 
pair transition density are confined in the nuclear volume, a proportional- 
ity relation B(P0) cx A 2 is expected^ in analogy with the B(E2) of the 
deformed rotor since the pair gap is a deformation parameter PEH We see 
in Fig. 5(b) that the proportionality B(P0) cx A 2 valid for 100 < A < 132 
is violated for A > 132 and especially A > 140, where the strength B(P0) 
significantly increases. The two-neutron transfer reaction on the neutron- 
rich isotopes (e.g. the Sn isotopes with A > 132) thus provides us a tool 
to probe the Cooper pair wave function in the low-density region far out- 
side the nuclear surface. It is predicted^ also that the isotopes 134_140 Sn 
with A = 134 — 140 exhibit a precursor phenomenon, i.e., an anomalously 
large two- neutron transfer strength of (t,p) type for the transitions to the 
excited + states (Fig. 5(b)). Recently two-neutron transfer experiment on 
the halo nucleus 11 Li has become available, and the crucial role of the pair 
correlation is demonstrated.^^ We wait for a future experiment using the 
radioactive beams of neutron-rich Sn isotopes with A > 132. 
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